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THEORY FOR COMPUTING SPAN LOADS AND STABILITY DERIVATIVES 
DUE TO SIDESLIP, YAWING, AND ROLLING FOR WINGS 
IN SUBSONIC COMPRESSIBLE FLOW* 

By M. J. Queijo 
Langley Research Center 

SUMMARY 

A method of computing span loads and some of the resulting aerodynamic derivatives 
for wings in sideslip, yawing, and rolling flight is derived. The method is applicable to 
wings of arbitrary planform, and accounts for compressibility effects. The span-load 
computations require that angle-of-attack span-load distributions be available for the 
wing under consideration. The theory developed is a consistent approach, based on the 
use of a vortex system, for determining the various wing parameters. 

INTRODUCTION 

The aerodynamic span loads and stability derivatives are important in relation to 
structural integrity and inherent stability of airplanes. The lifting surfaces (wing and 
tail surfaces) generally produce the predominant aerodynamic forces; therefore, much 
effort has been expended in developing methods of predicting the aerodynamics of lifting 
surfaces. Most of this effort, however, has been directed toward determination of aero- 
dynamic characteristics associated with angle of attack. Aerodynamics associated with 
other aircraft motions (rolling, yawing, pitching, and sideslipping) have been investigated 
to a lesser degree, and generally by somewhat cruder methods. This is particularly true 
for the low-speed regime. For supersonic speeds the nature of the governing equations 
is such that it has been possible to obtain equations for loads and aerodynamic derivatives 
for wings performing various modes of motion. (See refs. 1 to 4, for example.) In cer- 
tain limiting cases it is possible to use supersonic theory to predict subsonic character- 
istics. (See ref. 4, for example.) 

*A preliminary version of the material presented herein was included in a disserta- 
tion entitled ”A Theory and Method of Predicting the Stability Derivatives C^, C^, 

C n p, and Cyp for Wings of Arbitrary Planform in Subsonic Flow," offered in partial 

fulfillment of the requirements for the degree of Doctor of Philosophy in Engineering 
Mechanics, Virginia Polytechnic Institute, Blacksburg, Virginia, June 1963. 


There are a number of problems associated with attempting to predict aerodynamic 
characteristics of wings performing the possible modes of motion. One problem is that 
of finding an adequate mathematical model for the wing, and the second is that associated 
with solution of the equations which arise from use of the mathematical models. Various 
methods for predicting certain aerodynamic characteristics of unswept wings have been 
developed by a number of investigators, and numerous reports have been published from 
which certain characteristics can be obtained for specific wings. (See refs. 5 and 6, for 
example.) Three general approaches have been used in determining the aerodynamics of 
swept wings. These are: 


(a) Computations based on mathematical models associated with the use of vortices, 
doublets, or other concepts to represent the wing (refs. 7 to 11, for example). 

(b) Determination of approximate equations based on treating each wing semispan 
as one-half of an unswept wing. The fictitious "unswept" panels are skewed to simulate 
a swept wing (refs. 12 to 15, for example). 

(c) Development of design charts based on tests of a great number of wings with 
various sweep angles, aspect ratios, and taper ratios (ref. 16, for example). 


The first of these approaches is generally difficult and in some cases involves the 
solution of numerous simultaneous equations. For these reasons, only a few aerodynamic 
parameters ^primarily aerodynamic-center position, Cj^, and have been attacked 

by fairly rigorous methods. The use of high-speed computers to solve many simultaneous 
equations has permitted the numerical solution of equations better defining the wing bound- 
ary, but solutions generally have been obtained for angle-of-attack loading. 


The second approach has been quite successful in predicting trends, and with some 
modifications has been used to obtain good quantitative results for certain aerodynamic 
characteristics. (See ref. 15, for example.) 


The third approach is adequate for engineering data, provided the available data 
envelop the range of geometric variables of interest. Unfortunately, the amount of data 
available for some of the wing derivatives is very limited because of the scarcity of 
experimental facilities for determining such derivatives. 


The purpose of the present paper is to examine the problem of wing characteristics 
in subsonic compressible flow and to develop a consistent method for computing these 
characteristics. 


The method developed herein is a new approach for estimating span loads and the 
derivatives Cj , Cj , , Cn D , and Cy for wings of arbitrary planform in sub- 

p r p r p 

sonic compressible flow. It is based on a vortex representation of the wing which was 
first developed by the author for sideslipping wings in incompressible flow (ref. 17). 
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Results are generally applicable to the low angle-of-attack region, where the various wing 
characteristics vary linearly with angle of attack or lift coefficient. 

SYMBOLS 


aspect ratio, b 2 /s 
two-dimensional lift-curve slope 


wing span 


wing lift coefficient, 


lp V 2 S 





Cn 


pb 
9 — 
2V 

dCi_ 

9^ 

2V 

fC i 

9/3 


9^- 

2V 


Cy 


90 t 

wing lift-curve slope, 
rolling-moment coefficient, 


per radian 
M X 

|pv 2 sb 


M z 

yawing-moment coefficient, 

— pV 2 Sb 
2 


Fy 

side-force coefficient, 

1 P V 2 S 


3 


wing local chord 



8Cy 



2V 


c 

c 


c l 



wing average chord 

section lift coefficient, -- c - t f on 

|p y2c 

® C; 

section lift-curve slope, — - 

da 

three-dimensional section lift-curve slope for wing at angle of attack 
effective three-dimensional section lift-curve slope for rolling wing 

parameter for section lift, per unit lift coefficient 

parameter for incremental section lift due to sideslip, for a wing at angle of 
attack 

parameter for incremental section lift due to rolling 



c t 


D 

F Y 


parameter for incremental section lift due to yawing, for a wing at angle of 
attack 

wing root chord 
wing tip chord 

spanwise distance from wing root chord to center of rotation of yawing wing 
side force 
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h 

h 

h 

M 

m n 

M X 

M Z 

P 

r 
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V 

v N 

v 

X,Y 


x,y 


side force associated with a chordwise-bound vortex 
side force associated with quarter -chord-line vortex 
lift 

lift per unit length of chordwise-bound vortex 
lift per unit span of quarter-chord- line vortex 
lift per unit length of quarter -chord- line vortex 
free-stream Mach number 

Mach number of free stream normal to wing quarter-chord line 
rolling moment 
yawing moment 

rate of roll, radians per second 

yawing angular velocity, radians per second 

wing area 

wind velocity in x-direction 

free-stream wind velocity relative to wing center of gravity 
local velocity 

free-stream wind velocity normal to wing quarter-chord line 
wind velocity in y-direction 

longitudinal and spanwise reference axes, with origin at center of gravity 
distances along reference axes 
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*< I 


x ac chordwise distance between aerodynamic center (a.c.) and moment center 

(or center of gravity, c.g., in flight), positive when c.g. is upstream of a.c. 

x c x-distance to wing trailing edge 

x c /4 x-distance. to quarter-chord line 

spanwise position of centroid of the angle-of-attack span loading 

y radius of gyration of angle-of-attack span loading 

a angle of attack (or incidence), radians 

/3 sideslip angle, radians 

Ty circulation strength related to spanwise position 

circulation strength related to displacement along quarter-chord line 
6 local effective sideslip angle due to yawing, radians 

e infinitesimal displacement in spanwise direction 

A sweep angle of quarter -chord line, radians 

X wing taper ratio, c^jc r 

i ■ distance along the quarter-chord line, measured from wing root chord 

p mass density of air 

Subscripts: 

M compressible flow 

M=0 incompressible flow 

p part due to rolling 
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r 


part due to yawing 


a part due to angle of attack 

0 part due to sideslip 

A star (*) indicates that the quantity has been nondimensionalized by division 
by b/2; for example, y* = 


VORTEX SYSTEM 


The vortex system used in this analysis is for essentially a modified lifting-line- 
theory approach, and is illustrated in figure 1. The system consists of a bound vortex 
along the wing quarter-chord line and a bound-vortex sheet from the quarter-chord line 
to the wing trailing edge. Behind the wing trailing edge, the vortex sheet is in the direc- 
tion of the airstream. This system was applied to sideslipping wings in reference 17, 
and proved to be quite accurate in predicting the parameter Ci for a wide range of 

r 

wings in subsonic, incompressible flow. An appraisal of this vortex system indicated 
that it would be applicable to the estimation of other wing derivatives. 

The vortex system adopted for this study allows the possibility of lift generation by 
the bound vortices, which are: 


(a) The quarter-chord-line vortex, which extends across the entire wing span 

(b) The chordwise-bound vortices, which are parallel to the wing plane of symmetry 
and extend from the wing quarter-chord line to the wing trailing edge. 

The trailing vortex sheet behind the wing is made up of "free" vortices which are 
in the direction of the airstream and hence develop no lift. 


CIRCULATION DISTRIBUTION 

The strength of the chordwise-bound vortices is determined by the gradient of the 
strength distribution of the quarter-chord-line vortex; therefore, the lift distribution of 
the wing can be determined if the vortex strength distribution of the quarter-chord line is 
known and if the wind velocity components are known. The distribution of the wind veloc- 
ity components relative to the wing can be determined easily for each possible motion of 
the wing. The basic problem in determining wing load distribution, therefore, is that of 
determining the vortex (or circulation) distribution for all wing motions. The three types 
of motion considered in this paper are sideslip, yawing, and rolling. 
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Sideslipping Wing 


The sideslipping wing was studied in reference 17. It was argued therein that the 
circulation distribution in sideslip was essentially the same as that for a wing in zero 
sideslip. Thus zero-sideslip circulation distributions could be used for estimating aero- 
dynamic loads and the parameter Cj for wings in sideslip. The vortex system for the 

wing in sideslip and the direction of the airflow are shown in figure 2. 


Yawing Wing 

The bound-vortex pattern of the yawing wing is the basic pattern shown in figure 1. 
The trailing-free vortex sheet, however, is curved to match the airflow streamlines 
(fig. 3). An examination of the flow pattern over the wing (fig. 4) shows that there is a 
lateral velocity component resulting from the flow curvature, and that the magnitude of 
the component is a function of position on the wing. The wing therefore can be considered 
to be in sideslip, with the effective sideslip angle varying over the wing. The arguments 
presented with regard to circulation distribution of the wing in sideslip can be carried 
over to the yawing wing. An assumption of the present theory is that the circulation dis- 
tribution for a yawing wing is essentially the same as that for a nonyawing wing. 

Rolling Wing 

The problem of circulation distribution for the rolling wing is somewhat different 
from that for the sideslipping and yawing wing. In the case of the rolling wing, the local 
geometric angle of attack is increased by 



The primary cause of circulation, and hence the circulation itself, is altered by rolling. 
The net circulation of a rolling wing, therefore, is made up of: 

(a) Circulation due to the symmetric angle of attack 

(b) Circulation due to the antisymmetric angle-of-attack distribution associated with 

rolling velocity 

The assumed wing vortex system and the circulation distributions that have been 
discussed form the basis for the present theory for the computation of wing aerodynamic 
characteristics. 
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PRESENTATION OF RESULTS 


General Remarks 


General equations are derived in appendix A for the span loads and certain aerody- 
namic characteristics associated with sideslipping, yawing, and rolling for wings of arbi- 
trary planform. These equations are applicable in the low angle-of-attack region, where 
the characteristics vary linearly with angle of attack and lift coefficient. Span loads 
associated with sideslip, yawing, and rolling can be obtained by the methods presented 
herein only if the angle-of-attack span load is known. Such information is available for 
a wide range of wing geometry (ref. 18, for example). The angle-of-attack loading for 
odd- shaped wings can be obtained by application of the horseshoe-vortex method of ref- 
erences 19 to 22. 

The equations of this paper, for incompressible flow, are derived in appendix A. 
Compressibility effects are derived in appendix B. Some of the pertinent equations are 
summarized in appendix C. 

Many of the equations derived herein involve the spanwise position of the centroid 
and radius of gyration of the angle-of-attack loading of the wing semispan. The centroid 
position for angle-of-attack loadings has been determined for a wide range of wing plan- 
forms and is readily available in the literature. (See refs. 8, 18, and 21, for example.) 

The present theory appears to be the only one in the literature in which the radius of gyra- 
tion occurs as a factor in determining aerodynamic characteristics. It was necessary, 
therefore, to compute the radius of gyration y* for use in this paper. The values were 


and performing a mechanical integration to obtain y* by use of the following equation: 


A total of 160 such integrations were performed to obtain the desired coverage of wing 
planforms. Values of y* and of y* are plotted in figures 5 and 6 as functions of 
sweep, aspect ratio, and taper ratio for use in this paper and for general information. 
Data for the span loads were obtained from reference 18. 

Span-Load Distributions 

The span-load equations derived in this paper are applicable for the determination 
of the incremental load due to sideslip, yawing, or rolling velocity of a wing for which the 


determined by plotting the product 



against y* for a large number of wings 



( 1 ) 
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span load in symmetric flight is known. The span-load equations are summarized in 
appendix C. Equations (Cl) to (C3) are very general and should be used if the sweep 
angle varies across the wing span. In such instances it may not be possible to find in 
the literature the corresponding angle-of-attack loading of the wing. In these instances, 
as noted previously, the angle-of-attack loading can be determined in a straightforward 
manner (solution of simultaneous algebraic equations) by the method of reference 19, 20, 
or 21. 

If sweep is constant across the span, equations (C4) to (C6) are applicable and more 
convenient. In these cases the angle-of-attack span-load distribution can be obtained 
from a number of sources, such as reference 18. 

Most of the theoretical angle-of-attack span-load distributions available in the 
literature have been determined by use of the Weissinger method. It is of interest and 
significance that angle-of-attack span loadings determined by the use of horseshoe vor- 
tices (as in ref. 19 or 20) will closely approximate loadings from the Weissinger method 
if a minimum of about 20 horseshoe vortices are used to represent the wing. (See 
ref. 19.) The horseshoe-vortex method is applicable for representing wings having dis- 
continuities in sweep distribution (for example, M- or W-wings), where the applicability 
of the Weissinger method is doubtful. This horseshoe-vortex method was used in refer- 
ence 22 to determine the span loads of some M- and W-wings, and can be used for wings 
of arbitrary planform. The system of horseshoe vortices was also used in reference 17 
to obtain the loads on sideslipping wings. In such a case, the loads are obtained from 
either a mechanical integration or from summing terms of a series expansion (ref. 17). 


Aerodynamic Coefficients 

The equations derived for the various aerodynamic coefficients are presented in 
appendix C. Equations (C7) to (Cll) are very general and are applicable to wings of 
arbitrary planform. If sweep is constant across the semispan, and if the chord is a 
linear function of spanwise position, equations (C12) to (C16) are applicable and 
convenient. 


Equations (C12) to (C16) have been used to construct a series of figures from which 
numerical values of the various coefficients for a Mach number of zero can be obtained 
as a function of wing geometry. The effect of Mach number can be obtained through use 
of equations (C12) to (C16). The parameters for a Mach number of zero are presented 


in the following figures: in figure 7; C^C 

Cy^^CL in figure 10; and in figure 11. 


L in figure 8; C? in figure 9; 

r 
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DISCUSSION OF RESULTS 


Sideslipping Wing 

The span load due to sideslip can be computed by equation (Cl) or (C4). No other 
method appears to be available for estimating these loads. Experimental data are also 
almost nonexistent, except for the sideslip data given in reference 23 for a wing in incom- 
pressible flow. As shown in reference 17, these data compared well with the 
incompressible-flow equation of the present theory. 

The rolling-moment parameter associated with sideslip for a trapezoidal wing in 
compressible flow is given by equation (C12), repeated here for convenience: 


= 5 — + y.(tanA-«i^ 

Cl 2 A(1 + X) \ 


+ 0.05 


-Ay* 
2 y 


A 2 M 2 tan A 


( 2 ) 


H 1/2 


cos A 


- A 2 M 2 + 4 


2 + 


1/2 


cos A 


A 2 M 2 + 4 


The analysis of reference 17 shows that this equation is quite accurate for wings in 
incompressible flow. 

Very little theoretical work has been done to determine compressibility effects on 
the derivatives investigated in this paper. The Prandtl-Glauert transformation is appli- 
cable for symmetric flow, but there are a number of fundamental questions regarding its 
application if the flow about the wing is not symmetric. The Prandtl-Glauert transforma- 
tion can be applied either as a change in section lift-curve slope or as a change in wing 
geometry in symmetric flow. In unsymmetric flow, a geometric change in accordance 
with the Prandtl-Glauert transformation results in a distorted wing shape - that is, the 
transformed wing semispans are different from each other. This basic problem of appli- 
cation can be avoided, as was done in reference 14, by assuming that each wing semispan 
is one-half of a wing in symmetric flow. The compressibility corrections are then 
applied to each wing semispan and the loads of each semispan are used to obtain the wing 
characteristics. 

A somewhat simpler approach was used in reference 24, where the author simply 
chose to modify the section lift-curve slope as if the actual wing were in symmetric flow. 

It does not appear possible to obtain a direct evaluation of the estimated compres- 
sibility effects on the values of CiJ Cl for an isolated wing, since high-speed data are 
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generally obtained for wing-fuselage combinations. However, wing-fuselage data would 
be expected to have about the same trends as wing-alone data. Figure 12 shows experi- 
mental compressibility effects for several wing-body combinations and compares the 
results with those predicted by the wing-alone theory of reference 14 and by equation (2) 
of the present paper. The data and theories all show an increase in the magnitude of 
Cjc^ with increase in Mach number. The effect on wing-fuselage data is somewhat 

smaller than that predicted by the wing-alone theories. Reference 14 suggests an empir- 
ical method for accounting for the fuselage. However, as pointed out in that reference, 
additional effort is required to define more precisely the effects of the fuselage on 

%I C L- 

Yawing Wing 

The span loads due to yawing can be determined from equations (C2) and (C5). The 
author could find no other theory or experimental results to compare with these equations. 

The derivative obtained for the yawing wing is C j . For a trapezoidal wing in 

compressible flow it is given by 


i(l + tan2 A ) - » Ll± tan A + (y*) 2 tan A - 

2' ' 2Al+x 4^2\1 + \) ^ / \A1+X 2 )" > 


3 tan A 9(1 - X) 


2A(1 + X) a 2 (1 + x) 2 


y* + (MA . LLlAZj* + 3x ac + 


A l 


2A(1 + X) 4A 2 (1 + X) 2 


A 2 M 2 


-(y*) 2 tan 2 A + (y*^ Z (l + tan 2 A^ + x^ c y*tan A 


/^\2 


Q 

1/2 r 


1/2) 

( A ) - A 2 M 2 + 4 

\2 + 

2 

( A \ - A 2 M 2 + 4 

l 

Vcos A J 

l 

Vcos A I 

J 


( 3 ) 


Although this derivative is important in determining the dynamic lateral stability of air- 
craft, there is little published information on theoretical or experimental values of . 
Some incompressible-flow information is available from references 5 and 12. However, 
reference 5 is applicable only to unswept wings. Figure 13 shows comparisons of experi- 
mental data (incompressible flow) with the theory of reference 12 and with equation (3) of 
the present paper. Both theories fit the data equally well. 
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Reference 24 and equation (3) both predict an increase in with increase 

in Mach number, as shown in figure 14, but no experimental data are available for 
comparison. 


Rolling Wing 

The equation for additional load due to rolling velocity is, for a trapezoidal wing in 
compressible flow, 


r CCj> 

h. 


_l cc i 


ttA 


VcC 


L/M,a 


2 + 


1/2 


2 cos A 


A 2 M 2 


+ 4 


( 4 ) 


This type of loading is amenable to solution by the Weissinger method (ref. 25). Com- 
parisons of span load due to rolling as obtained by equation (4) and from reference 25 are 
shown in figure 15. The agreement between the two theories is quite good for incompres- 
sible flow. 

The span load distribution can be used to obtain the parameters , Cy , and 

r r 

C n , and these are given by the equations 

r 


C, = 




nA 


2 + 


2 cos A 


¥) 


1/2 


+ 4 


Cy 

2 = y*tan A 

c L 


(5) 


(6) 




tan 2 A 


+ x ac y*tan A' 


(7) 


Comparisons between equation (5) and experimental results for C? are shown in 

r 

figure 16. Generally there is good agreement between experiment and equation (5), par- 
ticularly at low aspect ratios. 
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The derivative Cy^ generally has no effect on the dynamic lateral stability of an 
airplane and therefore is not important in itself. However, its value is of importance in 
transforming certain of the other derivatives to various systems of axes. The basic ref- 
erence for theoretical values of this parameter in incompressible flow is reference 12. 
Theoretical and experimental results from reference 12 and values from equation (6) are 
shown in figure 17. Both theories generally underestimate this parameter. The present 
theory indicates no effect of compressibility on Cy^C]^ whereas reference 24 indicates 

that Cy p^L decreases with Mach number. There appear to be no available data to 

substantiate either theory. 

The wing contribution to C n p is important with regard to dynamic stability of 

airplanes. Several methods have been developed for predicting this derivative in incom- 
pressible flow. These methods are compared with each other and with experimental data 
in figure 18. Equation (7) indicates no compressibility effects on Cnpy^CL, whereas the 

theory of reference 24 indicates a decrease with increasing Mach number. However, 
there appear to be no experimental data available for determining compressibility effects 
on Cnp/CL* 


CONCLUDING REMARKS 

A theory and method are developed for computing span loads due to sideslip, yawing, 
and rolling for wings in subsonic compressible flow. The method is applicable to wings 
of any planform, provided the angle-of-attack load distribution is known. Such informa- 
tion is available in the literature for a wide variety of wing planforms, or can be obtained 
by straightforward methods. The span load (or circulation) distribution is used to deter- 
mine the stability derivatives due to sideslip, yawing, and rolling. Derivatives estimated 
by the method are compared with other theories and experimental results where available 
and applicable. It is not possible at this time to evaluate the present theory relative to 
others because of the limited amount of data available — particularly with respect to the 
derivatives associated with yawing and rolling. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 8, 1968, 

126-13-01-50-23. 
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APPENDIX A 


DERIVATION OF EQUATIONS FOR INCOMPRESSIBLE FLOW 


The vortex system used in this analysis includes a bound vortex along the quarter- 
chord line, a bound vortex sheet from the wing quarter-chord line to the wing trailing 
edge, and free vortices (along streamlines) behind the wing. Interaction of the relative 
wind with the bound vortices produces forces and moments. Equations for these forces 
and moments are derived in the following sections. In the derivations, reference is made 
to the right wing semi span unless otherwise noted. 


Sideslipping Wing 


The equations for span load and rolling moment for wings in sideslip were derived 
in reference 17. A few of the equations are included here for completeness and because 
they are required for the compressible-flow development. 


The incremental loading due to sideslip is given by 



(Al) 


and is antisymmetric over the wing. The resulting rolling-moment parameter is 



where the term 0.05 is a correction obtained in reference 17. 

For wings with straight leading and trailing edges the sweep angle is constant and 
c* can be expressed as a function of wing span: 


c* 


A(1 


7x)E " (1 " x)y *] 


(A3) 


Equation (Al) therefore can be written as 




tan A - 


S-4 

A(1 + A) l_ 


- (1 - A)y 


3 



dy* 


(A4) 
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and equation (A2) can be written as 


^/3 _ tan A 


c V cc z 

Jn 


'cc 7 


\cc L y a 


y*dy* + 2 f fl - (1 - x)y*|_X_Jd/" y*dy* + 0.05 (A5) 

, 2A(1 + X)JoL J dy* 


Equation (A5) can be integrated, with the result 


J1 I 

C L ' 2 


A(1 + X) 


+ y* tan A 


6 1 - X 
A 1 + X 


+ 0.05 


(A6) 


Yawing Wing 

The analysis of the yawing wing follows closely that of the wing in sideslip (given in 
ref. 17) but is more complex because the airflow velocity and direction relative to the 
vortices are functions of both the spanwise and the chordwise position on the wing. 

The total wind velocity at any point on the wing (see fig. 4) is given by 

V, = r(D - y) — i— - (A7) 

L cos 6 

The components parallel and normal to the wing plane of symmetry are 

u = -Vj cos 6 = -r(D - y) (A8) 


and 


v = -Vj sin 6 = -r(D - y)tan 6 = -rx (A9) 

For the right span of the wing, the component of velocity normal to the wing quarter-chord 
line is given by 


Vn = -u cos A - v sin A 
or, with substitution of equations (A8) and (A9), 

Vn = r(D - y)cos A + rx C/ / 4 sin A (A10) 
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\ 


If the velocity at the wing center of gravity is denoted by V, then 

V = rD 


and equation (A10) can be expressed as 

= (V - ry)cos A + rx c /4 sin A (All) 

The lift per unit length of the quarter-chord-line vortex is given by 

l ^ = pv N r^ (A 12) 

As discussed previously, a basic assumption of this analysis is that the circulation 
distribution is that due to angle-of-attack loading, and is not altered by yawing. The lift 
per unit length of the quarter -chord-line vortex therefore is given by 

h = ^* V N ( r £ ) a 

Similarly, the lift per unit span is given by 

h - <> v ( r y)a 

A unit length of quarter -chord-line vortex covers a span equal to cos A, so that 
l ^ = iy cos A. Using this equality, and noting that Vjq- = V cos A for the angle-of-attack 
case, it can be shown that 


(A13) 


(A14) 


{ r i)a ( r y)<* 

This simply indicates that in computing span loads the same circulation is used whether 
one deals with the free stream and the spanwise direction or with the quarter-chord line 
and the velocity normal to the quarter-chord line. Equation (A13) therefore can be 
written, per unit length of quarter-chord-line vortex: 

^ = PV N ( r y ) a 
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or, per unit span of the vortex: 

l y = p ^TA^y) a 

From equation (A14) it can be shown that 

( r y)« = \ Vc { c l)a 

Substituting equations (All) and (A16) into equation (A15) results in 


(A15) 


(A 16) 


h = i p £ v " ry) + rx c/4 



(A 17) 


per unit span. 

The lift of each chordwise-bound vortex must be found through an integration, since 
the lateral velocity varies along the length of the vortex. The lift is given by 


12 * C /4 pv ^2 d* 

x c dy 


or, with substitution of equations (A9) and (A16), 


fee. 


J 2 = -|prVcC L J Xc/4 : 


\cC 


La 


dy 


dx 


(A 18) 


The local section load coefficient is given by 


cc z ^ 1 ^y + ^2 

^ l p v 2 


\cC 


Lja,r --x, * pV ‘ 


or using equations (A17) and (A18) and normalizing with respect to the semispan gives 

_ f x c/4 rb x * 

IcCi 


cc 7 


cC 


L/Q!, r 


■ (‘ ■ § y * + w •“ A ) (srl - S'* § x * dx * (Al9) 


,cc l 

cCn 


dy 5 * 
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Since d 


r cc, 


'dy* is not a function of x, the integral of equation (A19) can be evalu- 


\cC L/d! 

ated easily, and equation (A19) becomes 


r cc. 


’ cc, 


cC 


L/Q!,r 


1 " !v( y * ' X c/4 tan A ) 


cc. 


icC 


L/a 


1 rb 

2 2V 


: c/4 


- ( x s) 


2l dl 


cC 


L/a 


dy 5 * 


The incremental load parameter due to yawing is given by 

d 

,\/ CC A 1 

4 c/4 


C ° l ^ = f-y* +x* /4 tan a)(^ 1 




Similarly, for the left wing semispan, 


cc i \ 

1 \ aC L/« 


cc L/a 2 d y* 


x c/4 


- ( x ?y 


r cc 7 


cc 7 


lcC.,4 


= (-y*-x; /4 tanA)(^- 


cC 


L/O! 


cc l/o! 2 


dy* 


2V/ 


( x S/4) 2 - ( x ?)‘ 


In general, the rolling-moment parameter for a wing is given by 

C i 


'l 
C L 


i r 1 / C V 


|y*dy* 


The rolling-moment parameter due to yawing is given by 

ly*dy* 


= .i r 1 

c l 4j -^ac L §/ 


(A20) 


(A21) 


(A22) 


(A23) 


(A24) 


Considerations of symmetry and the use of equations (A21) and (A22) with equation (A24) 
result in 


J - 1 Si' i ( y * - X c*/4 tan A ) (*£)„ + f( X J/4) 2 - H ) 


d( cc ' ' 

) 

\ C L/ 

la 

dy* 



J 


LZ“) y *dy* (A25) 
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Equation (A25) is perfectly general and somewhat complex for an arbitrary wing. If the 
wing leading and trailing edges are straight, however, some simplification is possible. 
In this case (see fig. 3), 


x S/4 = (y*-y*) tanA - x ac 


and 


x t = x c/4 " 4 c * 


Substituting these relations into equation (A25) results in 



(A26) 


Integration of equation (A26) is very lengthy and tedious. The solution of equa- 
tion (A26) is given by 


•i^l + tan^A^ - 


9 tan A 1 - A 
2A 1 + A 


27 | 

(\ - A 

4A 2 ' 

^1 + A 


~^\2 


(y*) 


— — — — tan A - — tan^A 
A 1 + A 2 


(y*)' 


3 tan A 

9(1 -A) 

V* 4- 

tan A 

3 1 - A 

-* , 3x ac 9 

2A(1 + A) 

A 2 (l + A) 2 

y + 

2 

A 1 + A 

“ 2A d - *> 4A 2 (1 + X) 2 


(A27) 


Rolling Wing 

In the analysis of the rolling wing the additional local angle of attack due to rolling 
must be considered. This additional local angle gives rise to an incremental circulation 
and therefore an incremental lift. Thus the net lift at any spanwise position is given by 
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This can also be expressed as 




(A28) 


(A29) 


where 


( c i ) and ( c l ) 

\ a) p \ a)a 


are the local three-dimensional section lift-curve slopes in 


rolling flight and at angle of attack, respectively. As a first approximation, it is assumed 
that the local lift-curve slope is equal to the lift-curve slope of the wing semispan. With 
this assumption, 


( c 0 « =Cl <* 1 

f'jp ” ( Cl o)p_ 


(A30) 


The additional lift due to rolling is antisymmetric over the wing; therefore, the effective 
lift-curve slope of the wing semispan in roll is approximately the lift-curve slope of a 
wing having the geometry (A and A) of the wing semispan. A convenient and accurate 
expression for wing lift-curve slope in incompressible flow is (from ref. 14) 


C L, 


Aa r 


a 


1/2 


(A31) 


2 + 


(COS A; 


+ 4 


Therefore, for the wing in roll the lift-curve slope is given approximately by 


( C Lq!)p 


TT a r 


1/2 


2 + 


1 2 cos A, 


+ 4 
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or, if a Q assumes its theoretical value of 2ir, 


( C]L “) P 


7rA 


2 + 


-, 1/2 


i2 cos A 


+ 4 


Using relations (A30), (A31), and (A32) with equation (A29) results in 

7TA 


< 


( c z)o-,p ( c z)q! 


2 + 


1 +• 


2 cos A 


+ 4 


1/2 


2V y 




(A3 2) 


(A3 3) 


Since the lift coefficient given by equation (A33) is all associated with the angle-of-attack 
type of loading, it can be directly related to circulation through equation (A16), so that 


ttA 


*\ 


r y = i Vc ( c i)« 


2 + 


1 + 


~ 1 1/2 


2 cos A 


+ 4 


£b 

2V 

J 


(A3 4) 


It is now necessary to consider forces and moments due to interaction of the veloc- 
ity components with the vortex system having the circulation distribution given by equa- 
tion (A34). The velocity components are the rectilinear velocity V and the velocity 
(normal to the wing plane) py. 

The interaction of the velocity V with the quarter-chord-line vortex results in a 
wing rolling moment given by 

r b/2 

M X = - ; pVr y dy (A3 5) 

-b/2 y 

or using equation (A34), nondimensionalizing, and taking the derivative with respect to 
pb/2V gives 
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C 



77 A 

2 


2 + 


From the definition of y* in equation (1), it is seen that 



(A3 6) 


Next the interaction of the rolling velocity py with the vortex system is consid- 
ered. This velocity is normal to the wing vortex plane. Since flow through the vortex 
sheet is not permitted, only the interaction of the additional velocity py with the 
quarter-chord-line vortex and the extreme wing-tip chordwise-bound vortices need be 
considered. 

The force caused by interaction of py with a chordwise-bound vortex is given by 


This force is in the wing (or vortex) plane and is in the y-direction (perpendicular to the 
chordwise-bound vortices). One approach to evaluating this force for the wing-tip vorti- 
ces is to perform the integration 



or, with substitution from equation (A16), 
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or 


n 


fee. 


\cC 


A(f Y )i=|ppVcC L 


h/a 


b/2 


yc 


f cc. 


u 


\cC 


h/a 


(b/2) -e 



This equation can be evaluated by parts, with the result 



1 =tPPVcC l < 


yc 


cc z \ 

CCl/q; 



O' 


b/2 



(b/2)-e 







Since all the functions involved in this equation are regular and well-behaved in the indi- 
cated limits, it is apparent that A^fy^ will approach zero as e approaches zero. 

Therefore, from the analysis made in this paper, the wing-tip vortices will contribute no 
resultant force or moment due to rolling velocity — that is, A^fyji = 0. 

The next possibility to consider is the interaction of the rolling-velocity components 
with the circulation on the quarter-chord line. The associated force is in the wing plane 
and is perpendicular to the quarter-chord line. The magnitude of the force is 

(fy) 2 = PPyFy (A3 7) 

per unit of vortex length. This force results in a side force which can be written as 


Fy = Vp 2v y * r y tanA 


(A3 8) 


per unit of span. A resulting side-force coefficient is given by 


pb/2 

c Y =y 


b/2 |pV 2 S 


dy 


(A3 9) 


Using equations (A34) and (A38) with (A39) and retaining up to first-order terms in 
pb/2V results in 


24 


APPENDIX A 


CY ^0 CL fe)a^ y * ta " Ady ’ 


The coefficient of side force due to rolling is found to be 

C Y r 


C L 


P _ 


= tan A 


I 1 


cc 


l 


0 V^L/ar 


y*dy’ t 


which integrates to 


■ D — 

-i- = y*tan A 


(A40) 


The force given by equation (A37) also results in a yawing moment which is given by 


r b/2 
-b/2 


M Z = J_ h/ ,( X c/4 Sin A - y C0S dy 


( J y) 2 


or, proceeding as before to form the yawing-moment coefficient, 


c r 
_±2 


£ 4 tan A ■ y *) : y * <iy * 


Integration of this equation results in 


5 = + 


C L 2 


(y*Y + (y *y - (y *) 2 


tan^A + x* c y*tan A ' 


(A41) 
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COMPRESSIBILITY EFFECTS 


The effects of compressibility on the span loads and derivatives can be estimated by 
proper treatment of the local span load. As a first approximation, the load coefficient at 
a particular spanwise position is assumed to increase in direct proportion to the total 
wing lift coefficient. At a constant angle of attack, therefore, 


( c l) M _ ( C l)m _ ( Cl o;)m 
( c l) M=0 ( C l)m= 0 ( C L Q! ) m=0 


Note that this assumption indicates no effect of compressibility on the center of pressure 
of the span load. 

A convenient equation for wing lift-curve slope, given in reference 14, is 


Aa r 


c l or 


1 - M 


N 


1/2 


7T 1 - M 


A 1/2 


N ) 


/ A ^ 

\ , 0 

/cos A J 

' 7T 2 I 

(l - M N 2) 


1/2 


(B2) 


Combining equations (Bl) and (B2) and letting a Q assume its theoretical value of 2 tt 
results in 



Equation (B3) can be used to modify the local span load and resulting derivatives to 
account for Mach number effects. 


Sideslipping Wing 

Reference 17 shows that for a sideslipping wing at an angle of attack in incompres- 
sible flow the local span-load parameter is given by 
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(1 + /3 tan A) - — c*j3 
a 4 



(B4) 


The Mach number perpendicular to the wing quarter-chord line is given by 

= M cos(A - (3) 

Combining equations (B3), (B4), and (B5) results in 


(B5) 


2 + 


cc 


cos A 


il/2 


+ 4 


c /M,j3,a 


1/2 


2 + 


cos A 


1 - M 2 cos 2 (A - /3)| +4 


-1 


cc 
c la 


cc l 


-i) (1 + p tan A) - -I c*/3 ^ e - ' a 

“ ' 4 r dy* 


(B6) 


The span-load increment due to sideslip is found by differentiating equation (B6) with respect to 0 and ignoring 
higher order terms in /3 (letting (3 approach zero). The result is 


CC A 

5(3 / M 


2 + 

/ a y 

\cos A ) 

2 1 
1 +4 

1/2 

O . 

( A ' 

\ 2 ( 

1/1 TV 

iT 2 rv /n n 9. A 1 i /I 

1/2 


icos A 


*1 tan A - - c* ^ C Ja 


cc 
c la 


dy’ 1 


A 2 M 2 (tan A)\ 2 + 

c a 


cos A 


+ 4 


1/2 


CO 

-a 


^ j ^1 - M^cos^a) + 4 


1/2 r 


r 


cos A 


M2cos2Aj 


+ 4 


1 / 2 ' 
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It should be noted that for zero sideslip, equation (B3) reduces to 


f V 

2 + 

/ A \ 2 + 4 

1/2 

( c l ) M,a _ 

T 

(cos A J 


( c l)o, 

2 + 

/ A 

(l - M 2 cos 

2 a) 

L-I T 

\cos A 

A / 


172 


(B 8 ) 


Using relation (B 8 ) with equation (B7) and dividing through by total lift coefficient results in 


cc. 


\^ C L7M \ cC L(/m,Q! 


/cc l\ , A 3 * d \P C L,/M 

' tan A - — c* 


dy* 


?J? l ) 

\ cC UM,a 


A 2 M 2 tan A 


cos A 


“1 1/2 


A 2 M 2 + 4 


2 + 


1_aJ ) 2 . a 2 M 2 + 4 
cos A/ 


1/2 


The rolling-moment parameter due to sideslip is (from ref. 17) 

'1 / cc. 


(B9) 


C 7 

Jl 

C T 


= .1C 

2 Jq 


7 


icC T 0/ 


y*dy* + 0.05 


or, using equation (B9), 
/C, 


Ife] ..if 1 

C L/M 2 *^0 


-*£< 


cc 


l 


VcC 


cc. 


L/M, 0 ! 


tan A - — c* 
4 


d& 

3 „* \ cC l/m,Q! 


dy* 


y*dy* + 0.05 


A 2 M 2 tan A 


f C L/M,a 


■ cos A 


A 2 M 2 + 4 


1/2 


2 + 


^— ) - A 2 M 2 + 4 


cos A 


1/2 


y*dy* 
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The correction factor 0.05, as discussed in reference 17, is independent of aspect ratio or sweep angle. Since com- 
pressibility effects can be related to geometric changes by the Prandtl-Glauert relations, and since geometry has no 
effect on the factor 0.05, it would appear that the factor is not a function of Mach number. 

Performing the integration as was done for the incompressible-flow case yields (for a trapezoidal wing) 


'L/M 


A(1 + X) 


+ y* tan A 


6 1 ~ X ' 
A 1 + \) 


+ 0.05 --y* 
2 J 


A 2 M 2 tan A 


Vcos A y 


A 2 M 2 + 4 


1/2 


2 + 


— — — ) - A 2 M 2 + 4 

^cos A/ 


1/2 


(Bll) 


Comparison of this equation with equation (A6) shows that 


U2 = Ufi 

: L/M \ C L/M=0 


1 -* 
" 2 Y * 


A 2 M 2 tan A 



1/2 " 

— 

1/2] 

/ \2 

/ \ ^ 

( 

( A J - A 2 M 2 + 4 

| 2 + 

( A ) - A 2 M 2 + 4 


\cos A/ 

L 

ycos A j 

J 


Yawing Wing 

The span load for a yawing wing in compressible flow can be determined in a manner analogous to that for the 
sideslip case. In yawing flight, however, the effective Mach number varies across the wing span. The velocity nor- 
mal to the wing quarter -chord line is given by equation (All), and the corresponding Mach number is 


M N = M 


cos A - rb 


|j^y*cos A - x*/ 4 sin h\ 


(B12) 
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Using equations (B3) and (B12) with equation (A20) results in 


2 + 


cc 


l 


\cos Ay 


+ 4 


1/2 


rb / * „* 

2V* 


i--^(y _ x c /4 tan A 


'cc 


l \ 1 rb 


lcC L l 2 2V 




x c/4 ) - 


cC 


L/ a 


dy* 


cC 


J 


L/ M,r,a 


2 + < —2 — m2 
cos A cos^A 


1/2 


cos A - -|^y*cos A - x^ 4 sin A 


+ 4 


(B13) 

Differentiating equation (B13) with respect to rb/2V, letting rb/2V approach zero, and proceeding as in the side- 
slip case, results in the following span-load parameter: 


'cc 


l 


cc, 


\ cC, 


rb 

2Vj 


= l-y* + xJ/ 4 tan A 


cc 


l 


icC 


M 


L/M,a 


c/4 


- (*c) : 


\cC 


L/M,^ 


dy* 


'cc 


l 


a2m 2 ^,y* + x*^ 4 tan Aj 


\cC 


L/M,0! 


cos A 


- A 2 M 2 + 4 


1/2 


2 + 


] - A 2 M 2 + 4 

cos A 


1/2 


(B14) 


The rolling-moment parameter due to yawing is found by substituting equation (B14) into equation (A24) and per- 
forming the integrations in the same manner as was done in the incompressible case. The result is found to be 
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'L/M 


M 1+tan2A ) • 


9 tan A 1 - X [ 27 / I - X ' 
2A 1 + X 4^2yl + Xy 


( y '*) 2 + ( lf 7 I tanA 4 tan2A )( r ) 


3 tan A 

9(1 - X) 

-* , ( tan A 3 1-1] 

2A(1 + X) 

A 2 (l + X) 2 

Y ^ 2 Al+xj 


, -* 3x ac 9 

x ;.y + xtt: — tt + 


4A 2 (1 + x) 2 


A 2 M 2 


~(y*) 2 t an ^A + ^y* ^1 + tan 2 A^ + x* c y*tan A 


cos A 


- A 2 M 2 + 4 


1/2 


2 + 


( - ) - A 2 M 2 + 4 

l cos A/ 


1/2 


Comparison of this equation with equation (A27) shows that 

,2 


'ly 


^r\ + A 2 M 2 


-Yy*^ tan 2 A + ^y* j ^1 + tan 2 A^ + x* c y*tan A 


ic l/m \ c l/m=o 2 


A ) - A 2 M 2 + 4 

mos A / 


1/2 


2 + 


) - A 2 M 2 + 4 


icos A I 


1/2 


Rolling Wing 


In compressible flow, equation (A28) becomes 


( c 0m,O!,p- ( c z)m,« 


! , fcKp 

( c ?)m,0! 


(B15) 


(B16) 
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The compressibility correction to p is based on the aspect ratio of a wing semi- 

span, as discussed previously; therefore, from equation (B2), 


( Cl o)m,p 


2 + 


7TA 


A '"'l-M N 2 + 4 


t2 cos A 


1/2 


and 


( c ;)m, p _ 


7rA — y* 
2V 


2 + 


2 cos A 


A Y, l- M n 2 J + 4 


1/2 


(B17) 


It can be shown that equation (B16) becomes 


flA 


2 + 


( c i) M,p,a - (^lj M,a 


1 +■ 


- A -) 1-M N 2 ) + 4 

2 cos A / ' N 


( c l) 


1/2 2V 


OT7- J 


M 




(B 18 ) 
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In rolling flow the velocity perpendicular to the wing quarter-chord line is given by 


V N 


2 

(V cos A) 2 + V 2 /^ (y*) 


1/2 


and the corresponding Mach number is 

Mn = M cos^A + 


-.1/2 


'§)V) : 


(B19) 


Substituting equation (B19) into (B18), differentiating with respect to pb/2V, and letting 
pb/2V approach zero results in the span-load parameter 


cc 7 


cc, 


—I 

2V/M 




2 + 


1/ A \ 2 A 2 M 2 


4 \cos A 


+ 4 


1/2 
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The circulation for the rolling wing in compressible flow is found by use of equa- 
tions (A16), (B18), and (B19), and is 


( r y) M,p,a--Z Vc { c l)M,a 


1 + 



7rA 

P b _ v * 

f( . \ 2 

/ , \2 

1/2 2V 

2 + |(o — ) 

^2 cos A/ 

1 - M 2 cos 2 A - M 2 ^tJ ^y*^ 2 

+ v 


( c l)m 



(B20) 


Substitution of equation (B20) into (A35), nondimensionalizing, taking the derivative with 
respect to pb/2V, and letting pb/2V approach zero results in the rolling-moment 
parameter 

c = -i rfe) - 

*p 2 Jq \cC l / m q, r . .o ~ ~ n 1 / 2 


2 + 


(y*) 2 dy* (B21) 


II A \ 2 A 2 M 2 


4 Vcos A 


+ 4 


Performing the integration results in 


C ; = 
Z P 


1 

2 


ttA 


2 + 


1/ A ) 2 A 2 M 2 


4 Vcos A 


+ 4 


1/2 


(y*)' 


(B22) 


Comparison of equations (A36) and (B22) shows that 


2 + 


(%)m- ( Cz p)m=o 


V2 cos Ay 


+ 4 


1/2 


2 + 


1/ A \ 2 A 2 M 2 


4 Vcos A 


+ 4 


1/2 


(B23) 


The side force due to rolling is found through use of equations (A38), (A39), and 
(B20), and is (for pb/2V approaching zero) 

Cy 

— — £ = y*tan A (B24) 
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Comparison of equations (B24) and (A40) indicates that there is no compressibility effect 
on the parameter Cy p ^C L . Similarly, it can be shown tlmt there is no compressibility 

effect on the parameter C n p /c L . 


34 


APPENDIX C 


Sideslipping wing: 


SUMMARY OF EQUATIONS 


General Equations for Additional Local Span Loading 
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Yawing wing: 
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Rolling wing: 
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Equations for Additional Local Span Loading for Swept Wings 


Sideslipping wing: 


cc, \ /cc 
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Yawing wing : 
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Rolling wing: 
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General Equations for Aerodynamic Derivatives 


Sideslipping wing : 
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Yawing wing : 


Rolling wing: 
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Equations for Aerodynamic Derivatives of Swept Wings 


Sideslipping wing : 
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Yawing wing : 
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Rolling wing: 
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Figure 2.- Complete vortex system for a wing in sideslip. 



Figure 3.- Complete vortex system for the yawing wing. 
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(a) x = o. 

Figure 5.- Spanwise location of centroid of angle-of-attack loading. 



47 













A, deg 



























-45 











































-30 















0 





























Aspect ratio, A 


(C) \ = 0.5. 
Figure 5.- Continued. 
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(d) X = 1.0. 
Figure 5.- Continued. 
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(e) x = 1.5. 
Figure 5.- Concluded. 
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(a) x = o. 

Figure 6.- Spanwise location of radius of gyration for angle-of-attack span load. 



(b) A = 0.25. 
Figure 6.- Continued. 
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(c) x = 0.5. 
Figure 6.- Continued. 
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A, deg 



(d) x = 1.0. 

Figure 6.- Continued. 
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(e) x = 1.5. 
Figure 6.- Concluded. 




Figure 7.- Variation of with aspect ratio, sweep, and taper ratio. M = 0; Xg C = 0. 
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Figure 7.- Continued. 



Figure 7.- Concluded. 




Figure 8.- Variation of Cj /c L with aspect ratio, sweep, and taper ratio M = 0; x| c = 0. 
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Figure 8.- Continued. 



Figure 8.- Concluded. 




Figure 9.- Variation of Cjp with aspect ratio, taper ratio, and sweep. 
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Figure 9.- Continued. 



Figure 9.- Concluded. 
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(b) X =0.25. 
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Figure 10.- Variation of Cy p jCL with aspect ratio, sweep, and taper ratio. M = 0; Xg C - 0. 
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(c) X = 0.50. 









Figure 10.- Continued. 



Figure 10.- Concluded. 
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Figure 11.- Variation of 


Cn p /C L 


with aspect ratio, sweep, and taper ratio. M = 0; 
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(c) X = 0.50. 















Figure 11.- Concluded. 
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Experiment, ref. 14 
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Figure 12.- 


Experimental and calculated effects of compressibility on 
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Eq. (3) 

Theory of ref. 12 



Figure 13.- Theoretical and experimental variations of C[ r with C L for several representative wings in incompressible flow. 
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Figure 14.- Estimated effect of Mach number on the parameter 



Figure 15.- Computed span-load distribution due to rolling velocity. M = 0. 
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Figure 16.- Calculated and experimental values of damping-in-roll parameter Cj p . M - 0; X - 1.0. 



Figure 17.- Theoretical and experimental variations of Cy p with Cl for several representative wings in incompressible flow. 
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Figure 18.- Theoretical and experimental variations of C n p with C|_ for several representative wings in incompressible flow. 
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